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Abstract 

The influence of external magnetic field h on a static conductivity of Mott- 
Hubbard material which is described by model with correlated hopping of electrons 
has been investigated. By means of canonical transformation the effective Hamil- 
tonian which takes into account strong intra-site Coulomb repulsion and correlated 
hopping is obtained. Using a variant of generalized Hartree-Fock approximation the 
single-electron Green function and quasiparticle energy spectrum of the model have 
been calculated. The static conductivity a has been calculated as a function of h, 
electron concentration n and temperature T. The correlated hopping is shown to 
cause the electron-hole asymmetry of transport properties of narrow band materials. 
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The achievements of the recent years in the field of strongly correlated electron sys- 
tems give us the opportunity to understand the properties of narrow-band materials, in 
particular those in which metal-insulator transition under the action of external influences 
(pressure, doping, temperature) is realized [1]. The strongly correlated electron systems 
demonstrate unusual transport properties [2]. For understanding of the physical mecha- 
nisms, which cause these peculiarities, the experimental and theoretical researches of the 
temperature dependence of conductivity are needed. The results concerning low-frequency 
behavior of conductivity are of the prior importance, because it gives the information 
about the scattering processes close to the Fermi surface. The theoretical investigation of 
conductivity a(u,T) are mainly concentrated in the limit T = 0. Behavior of the static 
conductivity cr(T) = a{oj = 0, T) at T > has not been studied sufficiently. 

Theoretical investigations of the optical conductivity in the Hubbard model [3] in the 
frameworks of the Kubo linear response theory [4] last for many decades, we note here 
the investigations by analytical methods: moment method [5], in composite operators 
method [6], in the mean-field theory [7, 8], in the perturbative theory method [9, 10] 
in the limit of weak interaction (|t| ^> U), in method of the memory function [11, 12], 
in the opposite limit {U ^> \t\). The conductivity has been intensively studied in one- 
dimensional Hubbard model [13]- [15], where the numerical results can be compared with 
exact ones obtained by Bethe ansatz application. 
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For the numerical investigation of the conductivity in Hubbard model the exact diago- 
nalization of finite clusters has been used [16] for 4 x 4 sites cluster, quantum Monte-Carlo 
method for 8, 10 sites [17], 3x3 sites [15], 8x8 sites [18], 12 x 12 sites [19]. The inves- 
tigations by numerical methods were carried out mainly for two-dimensional lattice (this 
is caused by the interest to the high-temperature superconductivity phenomena [20] in 
systems with CuO planes), the conductivity of three-dimensional system has been studied 
only in narrow interval close to half-filling at weak and intermediate interactions [21]. 

The optical conductivity has been studied also by the dynamical mean field theory 
(DMFT) in the limit of infinite spatial dimension [22]. Different DMFT equations solvers 
were used: iterated perturbation theory IPT [23], non-crossing approximation NCA [24], 
second order perturbation theory 20PT [25] . The authors have considered the symmetri- 
cal Hubbard model at half-filling close to metal-insulator transition (intra-atomic Coulomb 
repulsion U ~ 2w, where 2w is the band width) and have obtained a good agreement with 
experimental data for some Mott-Hubbard systems. But for realistic models the non-local 
(dependent on wave vector) contributions to self-energy and transport characteristics are 
important. Besides, investigations in the framework of DMFT consider mainly the case 
of half-band filling. 

In works [26]- [28] the essential importance of taking into account in the Hubbard 
model the correlated hopping of electrons and narrow-band model with non-equivalent 
Hubbard subbands has been proposed. In such a model the hopping integrals which 
describe translation movement of holes and doublons, differ one from the other and from 
the activation processes integral. Similar models have been studied intensively in recent 
years [29]-[31]. For the generalized model which take into account the non-equivalency of 
Hubbard subbands, reliable results for conductivity do not exist, so the analytical study 
of conductivity in the framework of the realistic models of the electronic subsystem is of 
importance. 

In this work we show that application of the variant of projection procedure for Green 
function allows to reproduce some peculiarities of static conductivity of narrow-band 
material in the limit of strong Coulomb correlation and investigate the effect of external 
influences like temperature change, doping, pressure and magnetic field. We apply our 
approach to the Hamiltonian, which, besides of intra-site Coulomb repulsion U, strong 
in comparison with inter-site hopping tij, additionally describes correlated hopping of 
electrons (influence of electron concentration n on the hopping processes) and show that 
it leads to the electron-hole asymmetry of conductivity and other characteristics. 

We write the Hamiltonian of correlated electron system in representation of Xf Hub- 
bard operators: 

H = H + H 1 + H[ + H ex , (1) 

where 

ff„ = E ( x i + x ") + U J2 X " + \nv ku 2 + pt B h £ Vs x°, 

is i ~ is 

H X = ^(n)XfXf + ^U^Xf'X?, 

ijs ijs 
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H 'i = £' ^VsXfXf + h.c.) . 



ijs 



H ex = -^m ((xt+xf) (.v; • .vj) . .vr.vr) . 



IJS 



Here operator Xj° l describes transition of site % from state |Z) to state \k), n is the chemical 
potential, U denotes the energy of intra-site Coulomb repulsion of electrons, J stands 
for the direct inter-site exchange interaction, k is the elastic constant, Vo is the initial 
volume of the crystal, N is the number of lattice sites, /is is Bohr magneton, h stands 
for the external magnetic field, n s = 1 if s =f and —1 otherwise. Translation processes 
are characterized by different hopping integrals, namely tij(n) = (1 + cra)(l — T\n)tij 
and tij(n) = (1 + au)(l — Tin — 2r 2 )tij are hopping parameters for holes and doublons, 
respectively; t'^in) = (1 + au)(l — Tin — T 2 )tij is parameter of the hopping of an electron 
between doublon and hole; correlated hopping parameters t 2 and T\ describe the influence 
of sites involved into the hopping process and neighbor sites, respectively; parameter a < 
takes into account the renormalization of bandwidth 2w = 1z\t i3 \ at strain u [32], z is the 
number of nearest neighbor to a site. 

We restrict ourselves to considering the strong correlation limit, namely U 3> w(n). 
In such a system at partial filling of the band the conductance is mainly due to electron 
hopping processes within the Hubbard subbands and interband hopping can be neglected. 
At these condition we apply the canonical transformation [28] to the Hamiltonian (1). 




(2) 



where 




(3) 



with L(ij) = t - 2 jj-^-- The operator S is taken to exclude the processes with pair hopping 
of holes and doublons in the first order in hopping parameter: 



H[ + [S; H ] = 0. 



(4) 



Finally, we obtain the effective Hamiltonian: 



H e ff — H + Hi + H ex + H ( 



(5) 



where 




(6) 



IJS 



with the indirect exchange interaction parameter J(ij) 



u 
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Using the projection procedure [33] for the case of n < 1 we obtain for the single 
particle energy spectrum: 

E s (k) = -/i - z.Jn a - zJn<j + a s t k (n) + (3 S , (7) 

where the correlated narrowing of the band and spin-dependent shift of subband center 
are 

2 — n + i] s m n 2 — m 2 

as = 2 + 2{2 -n + Vs my { } 

ft = - (2-n 2 + ,. ro) ?^< X "^ (9) 

respectively, and z is the number of nearest neighbors to a site. The respective results for 
n > 1: 

-Es(k) = — /i + C/ — zJn a — z.Jn 5 + agt^n) + /3 S , (10) 
_ _ n-M^m n 2 - m 2 

as " 2 + 2(n + ^m)' 1 j 

A = - 7 ^ , E*k(n)W 2 Xf) k , (12) 

describe upper Hubbard subband of halfbandwidth u)(n) = 

With use of the method of works [34, 35], we calculate the xx-component of static 
electronic conductivity a xx = a + a, where 



e 2 rz 



2Na *r. 

ijs 

is the conductivity of lower (0 — s)-subband, 



* = -^E<*W>*«(")> ( 14 ) 



is the conductivity of upper (|| — s)-subband. 

The magnetization can be calculated from the equation 



2h 



+ + zJ eff m\ _ sinh sinh 



1-n 



V © >/ sinh sinh (^^7) 

for n < 1 and corresponding equation in which the substitutions n — > 2 — n, w(n) — > w(n), 
« s — > <S S , /5 S — * /3 S are made. Here J e /j = J — J. Using these expressions we have 
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Figure 1: The concentration dependencies 
of the conductivity at correlated hopping 
t~i — r 2 = 0.1 in absence of the exter- 
nal magnetic field and exchange interac- 
tions: curve 1 corresponds to temperature 
e/w = 0.01, curve 2 - to Q/w = 0.1, 
curve 3 - to Q/w = 0.2. 



Figure 2: The concentration dependencies 
of the conductivity in the external mag- 
netic field at Q/w = 0.02, zJ eff /w = 0.02, 
t~i = T2 = 0: curve 1 corresponds to 
h/w = 0, curve 2 - to h/w = 0.005, curve 
3 - to h/w = 0.02. 



numerically calculated the static conductivity a xx as a function of electron concentration 
(Figs. 1, 2), temperature (Fig. 4), and magnetic field (Fig. 3). 

From Fig. 1 one can see that in the considered model with correlated hopping of elec- 
trons the conductivity provided by carriers in upper subband is lower than the conduc- 
tivity, provided by carriers from upper subband. This effect was discussed in work [35], 
it is a manifestation of the electron-hole asymmetry, inherent to real transition metal 
compounds. Other important feature is the change of current carrier type from metallic 
to semiconducting type in the vicinity of n = 2/3,4/3 and from semiconducting type 
to metallic one at n = 1. Increasing the correlated hopping we shift maxima of the 
conductivity closer to the half-filling. External magnetic field changes the concentration 
dependence of a qualitatively (see Fig. 2). The higher is the electron concentration, the 
less pronounced is the effect of applied magnetic field. At small concentration of electrons 
the band is fully polarized, in such a ferromagnetic system the conductivity is consider- 
ably lower than in paramagnetic state. If electron concentration increase, the decrease of 
magnetization leads to the increase of conductivity, a approaches its value in paramag- 
netic state. At the same time, due to the changes in the correlation band narrowing and 
the shift of subband center, the position of conductivity maximum changes from n = 2/3 
in saturated ferromagnetic state (Fig 2, 2) to n — 0.5 in ferromagnetic state (Fig 2, 3). 
The change of conductivity with increase of the applied field can be very sharp at low 
temperature, (Fig 3, solid curve). The increase of temperature leads to the decrease of 
a and its changes become more smooth. Plateau of a{h) dependence signifies, that at 
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low temperatures the fully polarized state is reached in very low field. Much higher field 
is needed to polarize spins at higher temperatures (Fig 7, long-dashed and short dashed 
curves). Sharp increases of a with increase of temperature is related to the transition 



Figure 3: The dependencies of the con- 
ductivity on magnetic field for the elec- 
tron concentration n = 0.3 and zJ e ff/w = 
0.01, T\ = t 2 = 0: solid curve corre- 
sponds to Q/w = 0.02, dashed curve - 
to Q/w = 0.04, dashed-dotted curve - to 

e/w = 0.06. 



Figure 4: The dependencies of the con- 
ductivity on temperature in the external 
magnetic field h/w = 0.01 at zJ e ff/w = 
and T\ = r 2 = 0: curve 1 corresponds to 
n = 0.3, curve 2 - to n — 0.35, curve 3 - 
to n = 0.45. 



from the polarized ferromagnetic state to paramagnetic one, in each state the conductiv- 
ity decreases with rise of temperature. Very small changes of electron concentration can 
greatly affect the temperature dependence of conductivity (see Fig. 4). 

In the framework of the considered model one can naturally introduce the notions 
of "wide" (lower) and "narrow" (upper) energy bands and, correspondingly, "light" and 
"heavy" current carriers with the effective masses 

m » = {-^) ' m * = (^k^J ' (16) 

where E s (\t) is the energy spectrum of current carriers in lower (s — 0) band, E s (k) is the 
energy spectrum of current carriers in upper (s— ||) subband. 

We have obtained that effective mass of heavy carriers can increase substantially 
with the increase of electron concentration. Because of the correlated hopping the 
heavy carriers can be realized in lower (s — 0) subband as well. It is worthwhile 
to note that the notion "effective mass of current carrier" has conditional sense, dif- 
ferent from the standard one, used in band theory. The definitions (16) are related 
to the expressions for band spectrum which describe the transitions between \s)- and 
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|0)-states as well as | fl) - an d |s)-states. Then m* and m* are effective masses 
of respective transitions; for the cases when subbands are almost empty or almost 
full, m* and m* can be interpreted as effective masses of electronic and hole states. 
In the paramagnetic state for the case n <C 
1 the value m* = ~ ~k\ can be 

identified as an effective mass of |s)-states, 
i.e. electrons, for the case n = 1 — e 
(e < 1) the value m* = 2a 2^ (n) | U=i is an 
effective mass of a hole. If n > 1, then 
for n — 1 + e we have translational mo- 
tion of doublons ("extra electrons") with 
masses m* = 2a2 ^ (n) | U=i, and for n = 2-e 
we have the effective mass of holes m* 2 = 
2a 2 |l(n)| \ n=2- From these results one can see 
that the effective mass of carrier in upper 
band can differ substantially from the ef- 
fective masses in the case when the con- 
ductivity is due to s — O-transitions. It is 

important to note that passing from the 

c j , . .j j i ,i Figure 5: The dependence of current car- 

regime of conductivity provided by the car- ~ , 

• i i j , ,i • i • riers effective mass m e ft/m on electron 
riers m lower band to the regime when it is JJ 1 u 

■j j u - -m + -j.- i-u a 4-- concentration: upper curves correspond to 
provided by s— \ ^-transitions, the effective , n 

, • , . , , Ti = r 2 = 0.3: middle curves - to Ti = 

mass increase stepwise at the point n = 1 ' 

(the behavior of the effective mass dependence "on , 't%e elec^ro^TEncen^r^ in 
Fig. 5). 

The effective mass dependence on the magnetization and external magnetic field is 
shown in Figs. 6 and 7, respectively. Rise of magnetization leads to the rise of difference 
in effective masses of spin- up and spin-down current carriers (Fig. 6). The slope of the 
dependencies m e ff(m) changes with the rise of m, too. This leads to the decrease of 
overall transport in ferromagnetic state, though effective masses of carrier with majority 
becomes lower. The results shown in Fig. 6 qualitatively agree with the corresponding 
plot of work [36], where Gutzwiller approximation has been used to calculate effective 
masses of current carriers. The correlated hopping, favoring localization, shift effective 
masses up. Different possible scenarios of m e ff(h) dependence are shown in Fig. 7. At 
high temperature m e ff changes monotonically (Fig. 7, dotted curve) while at low temper- 
ature system in unstable towards the transition to polarized state (Fig. 7, solid curve). 
Direct exchange interaction can stabilize ferromagnetically polarized state in the less then 
half-filled band even in weak magnetic field (Fig. 7, dash-dotted curve). In this paper we 
have used a model with correlated hopping of electrons to study the influence of external 
magnetic field, temperature and doping on a static conductivity of Mott-Hubbard mate- 
rial. In the regime of strong Coulomb interaction hybridization of the Hubbard subbands 
do not contribute essentially to the transport properties so we have excluded it from the 
effective Hamiltonian by means of canonical transformation. The single-electron Green 
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Figure 6: The dependence of current car- 
riers effective mass m e ///mo on the mag- 
netization at n — 0.2, Q/w = 0.02, 
zJ e ff/w = 0: solid curves correspond to 
T i — r 2 = 0; dashed curves - to T\ = r 2 = 



Figure 7: The dependence of current carri- 
ers effective mass m e ///rao on the applied 
magnetic field: solid curves correspond to 
zJ e ff/w — 0, Q/w — 0.02; dotted curves - 
to zJ e ff/w = 0, Q/w = 0.05; dash-dotted 
curves - to zJ e ff/w = 0.05, Q/w = 0.02. 



0.1. 



function and quasiparticle energy spectrum of the model have been calculated using a vari- 
ant of generalized Hartree-Fock approximation. This procedure has allowed us to obtain 
analytical expression for the band narrowing factor and relatively simple and transparent 
equations for calculation of spin-dependent shift of the band center, magnetization, static 
conductivity and effective mass of the current carrier. The results of our study generalize 
the results of works [6, 16], where the static conductivity of the Hubbard model has been 
calculated as a function of electron concentration n, for a wider class of systems, for which 
the correlated hopping should be taken into account. In the limiting case of the absence 
of correlated hopping our results agree with the concentrational dependencies obtained 
in the composite operator method [6], exact diagonalization [16] and Monte-Carlo simu- 
lations [18] results. The temperature dependence of conductivity in paramagnetic state, 
calculated in this work, agrees with the corresponding results of DMFT [24, 37]. In the 
magnetic field static conductivity reflects the changes of single electron energy spectrum 
through correlation narrowing of the band and shift of subband center. The temperature 
and concentration dependencies of o are governed by changes of system magnetization in 
external magnetic field. We have found that in the ground state the saturated ferromag- 
netic state is stable while at non-zero temperature magnetization has the concentration 
dependence which agrees with work [38]. Such behavior of magnetization leads to a{n) 
dependence with maxima at quarter and three-quarter fillings in distinction from param- 
agnetic ones, obtained in work [35]. At non-zero temperatures the sharp changes of cr(n) 
dependence are possible. It is due to complicated character of temperature dependence of 
band narrowing factor. Effective mass of quasi-particles appear to be spin-dependent and 
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substantially varies with magnetic field. These results are in agreement with the analysis 
of work [36] and experimental data of work [39] for heavy-fermion compounds. Taking 
into account the correlated hopping which is inherent to real narrow band materials al- 
lows us to describe electron-hole asymmetry of these processes which is observed in real 
materials. 
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